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Abstract 


This  report  presents  the  design,  development,  and  tests  of  a  high  frequency  thermoacoustic 
refrigerator.  It  was  developed  to  operate  at  5000  Hz  with  a  piezoelectric  driver  and  air  at  1 
atmosphere.  The  device  is  compact  even  though  it  was  operated  with  a  cavity  tuned  to  two  half¬ 
waves  of  the  drive  frequency.  The  stack  is  unusual  in  that  it  is  made  up  of  random  fibers  which  can 
provide  an  enormous  surface  area  for  interestion  with  the  sound  field.  At  each  end  of  the  stack  heat 
was  transferred  to  copper  heat  exchangers.  Substantial  cooling  was  observed  with  this  refrigerator. 
T emperature  difference  of  50°  C  were  produced  for  sound  levels  of  1 60  dB  at  the  stack.  The  cooling 
power  of  this  device  was  between  0.4  and  1 .0  watt.  Such  limits  are  attributed  to  a  thermal  boundary 
resistance  at  the  stack-heat  exchanger  interfaces  and  to  the  short  stack.  Experiments  showed  that 
limitations  in  the  stack  length  could  be  overcome  by  using  multiple  stacks;  this  approach  could  be 
extended  to  many  stages  of  refrigeration.  High  frequency  operation  leads  to  a  large  critical 
temperature  gradient  125°  C/cm  in  our  case,  which  opens  up  the  possibility  of  producing  very  low 
temperatures  in  a  suitable  geometry.  The  compactness  of  our  device  is  well-suited  for  high  pressure 
operation  and  for  rapid  cool-down  operation.  The  impressive  performance  of  our  device 
demonstrates  its  use  to  a  variety  of  refrigeration  applications  ranging  from  the  cooling  of  biological 
samples  to  cooling  of  components  in  high  speed  electronics.  Further  developments  of  this  device 
would  be  the  use  of  gas  mixtures  at  high  pressure  which  it  can  take  mechanically  due  to  its  small 
size. 


High  Frequency  Thermoacoustic  Refrigerator 


I.  INTRODUCTION 


Since  the  discovery  by  Merkli  and  Thomann  [1]  that  cooling  can  be  produced  by  the 
thermoacoustic  effect  in  a  resonance  tube,  research  has  concentrated  on  developing  the  effect  for 
practical  applications.  One  approach  was  to  increase  the  audio  pumping  rate.  Whilst  the  experiments 
of  Merkli  and  Thomann  used  frequencies  of  around  100  Hz,  Wheatley  and  collaborators  [2] 
successfully  raised  the  operating  frequency  to  around  500  Hz  and  achieved  impressive  cooling  rates 
in  their  refrigerator.  This  triggered  a  wide  range  of  very  promising  thermoacoustic  refrigerators. 
Here  results  will  be  presented  on  a  refrigerator  developed  for  operation  around  5,000  Hz.  The  report 
brings  out  the  design,  the  fabrication  and  tests  on  its  performance. 

Some  of  the  reasons  for  using  high  frequency  pumping  in  a  thermoacoustic  refrigerator  are: 

-  the  device  will  be  compact  as  its  dimensions  scale  with  the  wavelength  of  the  audio 
drive. 

-  quick  response  and  fast  equilibration  rate  for  small  devices. 

-  it  is  a  convenient  frequency  range  for  piezoelectric  drivers,  which  are  light,  efficient, 
and  cheap. 

-  some  parts,  such  as  heat  exchangers,  can  be  fabricated  using  photolithography  and 
their  film  technology. 

-  power  density  of  device  can  be  raised. 

-  many  applications  require  small  compact  refrigerators. 

An  important  element  in  the  operation  of  a  thermoacoustic  refrigerator  is  the  special  thermal 
interaction  of  the  sound  field  with  a  plate  or  a  series  of  plates  known  as  the  stack.  It  is  a  weak 
thermal  interaction  characterized  by  at  time  constant  given  by  (ot«1  where  <o  is  the  audio  pump 
frequency  and  t  is  the  thermal  relaxation  time  for  a  thin  layer  of  gas  to  interact  thermally  with  a  plate 
or  stack.  The  amount  of  gas  interacting  with  the  stack  is  determined  by  a  thermal  penetration  depth 
6k  given  by: 

A  =  (2k/co)1'2  (1) 

By  increasing  o>  the  condition  of  weakly  coupling  is  met  by  a  reduction  of  6k  and  hence  of  t.  Here 
k  is  the  thermal  difiiisivity.  The  work  of  pumping  heat  up  a  temperature  gradient  as  in  a  refrigerator 
is  performed  by  essentially  the  gas  within  the  penetration  depth.  The  amount  of  this  gas  will  have  an 
important  dependence  on  the  frequency  of  the  audio  drive. 

It  may  seen  that  conduction  of  heat  is  a  slow  process  and  hence  that  this  would  limit  the 
operation  of  a  high  frequency  refrigerator.  Such  objection  can  be  set  aside,  at  least  for  some  of  the 
frequencies,  because  heat  conduction  is  not  a  slow  process  when  small  distances  and  masses  are  in 
question  [3,4],  This  fact  will  be  applied  to  the  high  frequency  refrigerator  presented  here. 

II.  GOALS 


In  this  project  the  main  goals  were: 

-  to  develop  a  high  frequency  thermoacoustic  refrigerator.  The  high  frequency  refers 
to  a  pumping  rate  of  5  kHz  which  was  chosen  arbitrarily.  It  is  about  a  factor  of  10 
above  the  frequencies  used  by  Wheatley  and  collaborators;  it  falls  within  a  frequency 


range  where  relatively  cheap  hi-fi  equipment  can  be  used,  and  it  falls  within  a 
frequency  range  where  piezoelectric  drivers  are  efficient. 

-  to  develop  a  device  which  can  be  adapted  for  miniaturization. 

-  to  study  the  performance  of  a  small  refrigerator  and  the  features  peculiar  to  its 
compactness. 

-  to  develop  a  simple,  compact,  and  cheap  device  which  can  be  used  for  cooling  small 
electronic  components  and  small  biological  systems. 

-  to  explore  the  field  of  thermoacoustics  at  high  frequencies  with  emphasis  on  new 
concepts  and  new  materials. 

Ill  APPROACH 


Because  the  frequency  range  for  our  refrigerator  is  an  order  of  magnitude  higher  than  in 
conventional  thermoacoustic  refrigerators,  the  approach  and  some  of  the  components  are  different 
from  the  500  Hz  range  refrigerators.  In  fact,  in  some  of  the  features,  high  frequency  operation  leads 
to  a  simplification. 

As  a  starting  point  the  refrigerator  was  developed  with  the  following  main  features: 

-  working  fluid:  air  at  1  atmosphere. 

-  driver:  piezoelectric  device. 

-  resonator:  cylindrical. 

-  stack:  random  fibers 

-  heat  exchangers:  copper 

This  approach  was  chosen  for  its  simplicity  realizing  that  it  is  far  from  ideal.  Its  is  a  first 
approximation.  As  the  development  progressed  changes  were  made  to  improve  the  performance. 
An  analysis  will  now  be  presented  of  each  component  so  as  to  show  the  advantages  of  our  approach 
as  well  as  some  of  the  limitations.  Fig.  1  shows  the  basic  set-up. 

Piezo  Driver 


Stack 

Fig.  1  Basic  high-frequency  refrigerator. 

A.  ACOUSTIC  DRIVER 


We  selected  a  piezoelectric  driver  of  the  bimorph  type,  an  example  of  one  being  the 
Motorola  KSN  1046,  horn-loaded  for  better  impedance  matching.  This  model  was  chose  mainly 
because  of  its  high  sensitivity  and  broad  frequency  response.  Some  of  its  characteristics  are: 

-  mass:  1.3  g 


-  sensitivity:  ~  95  dB/watt/m  but  it  can  vary  by  a  few  decibels,  depending  on  the 
unit. 

-  frequency  response:  4-27  kHz. 

Fig.  2  shows  the  frequency  response  as  provided  by  the  manufacturer:  there  are  large  variations  with 
each  unit.  The  horn  cone  has  a  maximum  diameter  of  4  cm.  This  driver  efficiency  can  be  as  high  as 
70-90%,  depending  on  the  load. 
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Fig.  2.  Typical  frequency  response  of  bimorph  driver. 

This  type  of  driver  has  ideal  characteristics  for  use  in  a  high  frequency  refrigerator. 
Dissipation  power  losses  are  very  small  since  a  piezoelectric  is  a  capacitor  (ours  has  a  capacitance 
C  of  145  nano  Farads)  whose  losses  come  from  the  hysteresis  behavior  of  the  dielectric.  Compared 
to  an  electromagnetic  driver  whose  voice  coil  has  typically  ~  8  ohms  resistance,  its  series  equivalent 
resistance  representing  losses  is: 

R„  =  (6.28  fC  tan  6)'1 

being  of  order  1-10  kilo-ohms,  depending  on  the  material.  Here  tan  6  is  the  dissipation  factor.  Since 
dissipation  power  is  V2/^  it  is  much  smaller  for  the  piezo  driver  than  for  the  regular  electromangetic 
driver.  The  piezoelectric  driver  is  a  voltage  device  while  an  electromagnetic  driver  is  a  current 
device. 


In  a  bimorph,  two  piezoelectric  discs  are  bonded  together  on  each  side  of  a  brass  shim;  they 
change  lengths  in  opposite  direction  with  applied  voltage  and  this  causes  a  large  bending  action  which 
is  coupled  to  a  cone  diaphragm.  This  device  behaves  like  a  bimetallic  strip  which  flexes  upon  heating. 
Fig.  3  shows  a  bimorph  driver. 


Fig.  3.  Bimorph  Driver. 

Being  very  light,  this  device  is  attractive  for  the  present  applications.  One  of  its  limitations  is  the 
maximum  electrical  power  input  that  it  can  handle,  being  ~  5  watts  before  causing  mechanical  damage 
to  the  bimorph. 

B.  STACK 

At  first  a  conventional  stack  was  tried;  it  consisted  of  a  narrow  strip  of  Mylar  sheet  to 
which  nylon  fishing  line  spacers  were  glued  and  then  it  was  rolled  into  a  cylindrical  shape,  as 
developed  and  described  by  T.  Holler  [5],  This  was  unsuccessful  here  because  for  a  stack  length  Ax 
of  0.5  cm  it  was  hard  to  maintain  uniform  spacing  but  also  it  was  difficult  to  make  good  thermal 
contact  to  a  heat  exchanger  at  each  end  of  the  stack.  This  was  abandoned  for  a  novel  approach, 
which  did  not  have  the  above  problems. 

We  chose  a  random  element  stack  using  cotton  wool.  This  is  shown  in  Fig.  4.  It  was 
pressed  so  as  to  have  a  length  Ax  =  0.5  cm.  Some  of  characteristics  of  this  stack  are: 

-  density  of  cotton  wool:  0.08  g/cm3 

-  thermal  conductivity  of  each  fiber:  0.04  W/m  °C 

-  average  diameter  of  each  fiber:  12  pm 

-  number  of  fibers  in  stack,  3  cm  in  diameter:  4x  106 

-  effective  total  perimeter  II  of  fibers  ~  1 5 1  m 

-  effective  cross-sectional  area  for  heat  pumping,  II6k:  7.5 x  1 0"2  m2 

-  total  active  area  of  stack  exposed  to  sound  field:  7.5*  103  cm2. 

The  above  quantities  are  maximum  values,  as  in  some  cases,  randomness  of  the  fibers  was  not  taken 
into  account. 
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Fig.  4.  Random  Fiber  Stack. 

C.  HEAT  EXCHANGERS 

Two  types  of  heat  exchangers  were  used,  both  were  made  of  copper. 

-  a  photolithographically  prepared  heat  exchanger  with  square  holes,  0.5  mm  x  0.5 
mm  and  solid  to  solid  spacers,  0.8  mm  x  0.8  mm.  This  exchanger  provided  a 
sound  transparency  of  25%.  Its  diameter  was  3.4  cm  and  it  was  fabricated  from 
copper  sheet  0.3  mm  thick.  Fig.  5  shows  the  heat  exchanger.  The  cold  and  the 
hot  heat  exchangers  had  the  same  geometry. 


Fig.  5.  Photolithographically  prepared  copper  heat  exchanger. 

-  a  copper  screen,  flattened  by  a  press,  with  square  holes  0.8  mm  x  0.8  mm  and  1.2 
mm  for  adjacent  wire  to  wire  distance.  Its  diameter  was  3.4  cm  and  it  had  a 
thickness  of  0.5  mm.  Sound  transparency  was  44%.  To  improve  heat  transfer  at 
the  hot  heat  exchanger  (since  it  handles  more  heat  than  the  cold  one)  it  was 
thermally  anchored  to  a  large,  0.5  cm  thick,  copper  heat  exchanger. 


Although  thin,  the  heat  exchangers  were  eveloped  to  maintain  heat  flows  of  ~2  watts 
without  creating  a  substantial  AT  across  the  heat  exchanger  (AT  is  less  than  0. 1  °C).  The  mass  of 
each  heat  exchanger  was  1.2  g.  Mechanical  contact  was  made  between  the  heat  exchangers  and  the 
cotton  wool  stack  and  this  provided  the  thermal  contact. 

D.  WORKING  FLTJTD 

As  mentioned  earlier,  we  started  with  air  at  1  atmosphere  as  the  working  fluid,  realizing 
all  along  that  this  is  not  the  optional  working  fluid.  Important  characteristics  of  air  are: 

-  thermal  conductivity:  0.26  mW/cm/°C 

-  density  at  1  atmosphere;  20°C:  0.00121  g/cm3 

-  viscosity  at  20°C:  18.1  ppoise 

-  speed  of  sound  at  20°C:  344  m/sec 

-  thermal  penetration  depth  at  5  kHz:  0.05  mm 

-  viscous  penetration  depth  at  5  kHz:  0.035  mm 

-  Prandtl  number:  0.707 

Better  performance  is  expected  in  a  mixture  of  Argon  and  Helium. 

For  the  specific  mixture  Ar0  36He0 M,  some  of  the  characteristics  which  show  more  promise  for  the 
refrigerator  are: 

-  thermal  conductivity:  0.09  W/m/K 

-  Prandtl  number:  0.351 

-  speed  of  sound  at  20  °C:  497  m/s 

E.  RESONATOR 

A  simple  geometry  was  chosen  for  the  resonator,  that  of  a  cylinder  closed  at  both  ends, 
with  a  drive  at  one  end.  Such  tube  resonator  was  used  in  2  different  ways,  after  modifications: 

(1)  as  a  1/2  wave  resonator  tuned  to  5000  Hz. 

(2)  as  a  double  half-wave  resonator  tuned  to  5000  Hz  (i.e.  the  half-wave  part  was 
tuned  to  5000  Hz  and  the  resonator  contained  one  full  wave).  This  offers 
interesting  advantages  (one  being  a  higher  Q)  which  will  be  discussed  in  the  next 
sections. 

Apart  from  the  stack  and  heat  exchangers,  acoustic  dissipation  in  the  resonator  is  due 
mainly  to  viscosity  and  heat  conduction  losses  on  the  surfaces  of  the  resonator.  In  addition  nonlinear 
dissipation  caused  by  circulation  effects  and  turbulence  will  come  in  at  high  sound  intensities. 
Viscosity  dissipation  can  be  represented  by  a  surface  resistance  R,  which  is  given  by: 


Rs  =  ‘/*(2np  U>t  (3) 

which  at  room  temperature  for  air  becomes  R,  ~  0.008  f\  Here  p  is  the  coefficient  of  friction.  Heat 
conduction  losses,  also  for  air  at  1  atmosphere  and  room  temperature  can  be  written  as: 

H,  -  (lO-'Vpctf'4  (4) 


Such  losses  will  have  an  immediate  impact  on  the  quality  factor  of  the  resonator  [6]  defined  as: 


Q  =  O)om/pc0 


(5) 


IV  ACHIEVEMENTS 


This  section  consists  of  3  parts:  experiments  and  tests,  analysis  of  performance,  and  discussion 
of  potential  of  approach. 

A.  EXPERIMENTS  AND  TESTS 


The  high-frequency  refrigerator  was  developed  to  test  the  concepts  of  thermoacoutics 
applied  here,  to  use  the  results  as  guideline  for  further  improvements,  and  as  a  test  unit  for  further 
developments  in  high  frequency  refrigeration.  Basic  tests  of  our  device  consisted  of  tuning  the  driver, 
adjusting  the  sound  intensity  level,  and  following  the  cooling  with  a  thermocouple  across  the  stack. 
Although  the  piezoelectric  bimorph  drivers  are  broad-band  in  frequency,  a  certain  amount  of  tuning 
was  necessary  to  take  advantage  of  the  small  resonances  with  their  increased  sound  intensity  level. 
The  operating  frequency  was  between  4  and  5  kHz;  the  corresponding  wavelength  in  air  at  1 
atmosphere  ranged  from  8  to  6.8  cm. 

Hence  a  V2-wave  resonator  at  5,000  Hz  would  be  ~3.4  cm  long.  This  type  of  resonator 
provides  the  opportunity  to  make  a  compact  refrigerator.  We  have  concentrated  our  efforts  mainly 
on  the  second  type  of  resonator,  the  double  half-wave  resonator.  It  is  tuned  to  -5000  Hz  but  it  is 
twice  as  long  as  the  Vi  wave  resonator  since  it  contains  2  half-waves  of  the  same  wavelength  as  the 
Vi  wave  resonator.  This  is  shown  in  Fig.  6  with  the  stack  and  heat  exchangers  positioned  at  the 
appropriate  position  with  respect  to  the  pressure  standing  wave  in  the  resonator. 


Piezo  Driver 
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Fig.  6.  Half-wave  and  double  half-wave  resonators. 

The  geometry  of  the  double  Vi  wave  resonator  provides  the  option  of  having  2  or  more  stacks  which 
could  be  eventually  connected  in  tandem.  Tests  on  such  arrangement  will  be  discussed  later. 

(i)  Tests  with  Vi  Wave  Resonator. 

As  the  temperature  difference  AT  across  the  stack  is  a  function  of  the  position 


of  the  stack  in  the  acoustic  standing  wave,  a  positioning  device  was  incorporated 
in  the  device.  This  is  shown  in  Fig.  7.  It  allowed  us  to  study  AT  across  the  stack 

Piezo 


Fig.  7.  Device  for  positioning  the  stack  and  heat  exchangers. 

as  a  function  of  its  position.  At  some  point,  the  temperature  change  due  to  the 
pressure  change  of  the  sound  field  is  balanced  out  by  the  fluid  displacement  in  a 
temperature  gradient  and  this  leads  to  a  critical  temperature  gradient  Afcri,  It  is 
defined  as  [7]: 


4r" '  TF  T  “  ("1) 


where  y  is  the  ratio  of  isobaric  to  isochoric  specific  heats,  Tm  is  the  mean 
temperature  of  the  fluid,  X  is  the  radian  length,  p  is  the  thermal  expansion 
coefficient,  and  x  is  the  stack  position  relative  to  the  pressure  antinode.  We  have 
investigated  the  importance  of  the  position  of  the  stack  relative  to  the  acoustic 
standing  wave;  this  is  shown  in  Fig.  8,  with  the  spatial  dependence  normalized  to 
the  sound  radian  length.  30 
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Fig.  8.  Temperature  difference  across  stack  as  a  function  of  stack 
position  in  standing  wave. 


The  results  show  how  critical  the  positioning  of  the  stack  is  and  how  the  direction  of 
the  pressure  gradient  in  the  acoustic  standing  wave  determines  the  sign  and  magnitude  of  AT.  In  this 
case  the  stack  length  Ax  was  0.5  cm. 

Having  established  the  position  of  maximum  AT  the  stack  was  fixed  at  that  position 
so  that  the  cooling  could  be  studied  as  a  function  of  sound  power  and  time.  Using  calibrated 
miniature  electret  microphones  in  the  body  of  the  refrigerator  the  sound  level  was  determined  at  the 
position  of  the  stack.  Atypical  sound  intensity  was  156  dB  and  this  corresponds  to  0.4  W/cm2.  For 
a  3  cm  diameter  stack  this  gives  an  input  acoustic  power  level  of  ~2.5  watts.  At  maximum  power 
from  the  driver,  the  temperature  difference  AT  between  the  hot  and  the  cold  end  of  the  stack  reached 
typically  39°C.  This  would  not  last  very  long  as  the  driver  invariably  cracked  due  to  the  large  stress 
on  the  bimorph  at  such  high  intensity  levels.  Fig.  9  shows  a  typical  cooling  curve.  Temperature 
differences  were  measured  by  a  copper-constantan  thermocouple. 


t  (sec) 

Fig.  9.  Temperature  difference  across  stack  in  the  presence  of 
intense  sound  at  ~5000  Hz. 

When  the  sound  power  switched  off,  the  temperature  difference  AT  across  the  stack  decayed  to  zero 
with  a  characteristic  time  constant.  The  magnitude  of  the  AT  achieved  in  each  case  depended  on  the 
quality  of  the  driver  and  how  well  it  was  tuned  for  highest  sound  intensity.  In  a  batch  of  drivers,  the 
sound  level  sensitivity  could  vary  by  up  to  10  dB. 

(ii)  Tests  with  double  Vj  wave  resonator  and  2  stacks. 

Here  we  incorporated  2  stack-heat-exchanger  units  at  appropriate  positions 
in  the  double  Vi  wave  resonator  refrigerator.  A  thermocouple  monitored  the  AT 
across  each  stack-heat-exchanger  unit  as  the  sound  level  was  raised  and  then 
maintained  at  the  maximum  level.  At  160  dB  sound  level,  one  stack  produced  a 
AT  of  25°C  while  the  other  one  produced  35°C  at  the  same  time.  The  difference 
is  probably  due  to  inaccurate  positioning  of  one  of  the  stack-exchanger  units  in 
the  resonator.  It  is  important  to  note  that  each  unit  was  thermally  isolated  from 
the  neighboring  one.  The  results  are  interesting  because  the  two  units  could  be 
attached  thermally  in  tandem  for  further  overall  temperature  reduction. 


(iii)  Tests  with  the  Double  V2  Wave  Resonator  and  1  Stack 

In  this  series  of  test  only  1  stack-heat-exchanger  unit  was  used.  When  tuned 
as  to  stack  position  and  maximum  sound  power  level,  a  AT  of  50°  C  could  be 
reached,  typically.  The  stack  was  located  just  before  the  last  pressure  antinode 
away  from  the  driver.  A  better  performance  was  achieved  due  to  the  higher  Q  of 
the  resonator  (compared  to  the  single  V2  wave  resonator)  as  a  substantially  larger 
AT  was  achieved.  As  mentioned  before,  the  limiting  factor  in  the  achieved  AT 
was  the  performance  of  the  driver.  We  used  the  drivers  at  their  limit. 

(iv)  Tests  of  Double  V2  wave  resonator  and  1  stack  and  a  Load. 

With  the  encouraging  results  obtained  above,  the  thermoacoustic  device  that 
we  have  studied  was  incorporated  into  a  refrigerator  by  anchoring  the  hot  heat 
exchanger  to  room  temperature  and  letting  the  cold  heat  exchanger  drop  in 
temperature.  The  standard  heat  exchanger  at  the  hot  side  was  thermally  anchored 
to  a  large  ribbed  copper  heat  exchanger  which  was  held  at  room  temperature. 
The  portion  of  the  refrigerator  from  the  hot  exchanger  down  was  surrounded  by 
a  vacuum  jacket  for  thermal  isolation.  This  allowed  the  cold  heat  exchanger  to 
cool  the  portion  of  the  resonator  in  contact  with  it  (20  cm3  volume);  however  a 
AT  of  50  °C  could  not  be  achieved.  Instead  a  AT  of  25  °C  was  reached  and  this 
deteriorated  with  time.  A  thermal  bottleneck  was  encountered  at  the  stack-hot 
heat  exchanger  interface  and  the  temperature  of  the  stack  at  the  hot  and  rose 
above  room  temperature.  Low  heat  transfer  at  interfaces  limited  the  performance 
of  our  refrigerator. 

(v)  Refrigerator  using  gas  mixtures  and  high  pressure  operation. 

The  experiments  and  tests  with  air  as  the  working  fluid  provided  us  with  basic 
information  on  a  high  frequency  thermoacoustic  refrigerator  and  how  it  could  be 
improved.  One  direction  still  not  addressed  in  this  project  was  the  optimal  choice 
of  working  fluid.  Since  a  limiting  factor  is  the  viscous  boundary  larger 
characterized  by  a  viscous  penetration  depth  8V  it  is  important  to  chose  a  fluid 
where  this  depth  is  less  than  the  thermal  penetration  5k,  i.e.  a  fluid  with  a  low 
Prandtl  number.  Hence  we  chose  a  mixture  of  64%  He  and  36%  Ar  whose 
Prandtl  number  is  0.3507  and  where  the  speed  of  sound  is  497  m/sec  [8], 
Compared  to  air  this  required  a  scaling  factor  of  1 .4  in  size  to  keep  the  resonance 
at  the  same  frequency  as  for  air.  Tests  were  made  with  this  mixture  at  1 
atmosphere.  Preliminary  results  were  encouraging  as  a  AT  of  30°C  across  the 
stack  was  easily  achieved  for  a  stack  of  fixed  position  (and  not  necessarily  located 
at  optimal  position). 

All  our  tests  were  performed  with  working  fluids  at  1  atmosphere. 
Improved  performance  is  expected  when  the  fluid  is  at  higher  pressures  because  of  scaling  similitude 
principles  [9]  and  because  this  provides  superior  impendance  matching  between  the  driver  and  the 
fluid.  In  fact,  we  are  in  an  excellent  situation  for  working  at  high  pressure,  since  a  small  refrigerator 
is  structurally  strong  enough  to  withstand  very  high  pressures.  Such  approach  remained  in  the 
planning  stages. 


B.  ANALYSIS  OF  PERFORMANCE 

(i)  Maximum  AT 

The  maximum  temperature  difference  that  can  be  produced  across  a  stack  results  from 
a  competition  between  the  temperature  change  due  to  an  adiabatic  pressure  change  of  the  working 
fluid  and  its  displacement  along  the  stack  which  has  a  temperature  gradient.  When  the  temperature 
rise  due  to  an  adiabatic  compression  is  greater  than  the  temperature  rise  due  to  the  displacement 
along  a  temperature  gradient  of  the  stack,  the  engine  works  as  a  heat  jump.  Conversely,  the  engine 
works  as  a  prime  mover.  The  critical  gradient  AT^  given  by  Eg.  6  separates  the  2  regimes.  For  our 
case  ATcrit=120°C/cm  and  this  determines  the  maximum  AT  which  can  be  achieved.  For  a  stack  of 
length  Ax=0.5,  cm  as  in  our  case,  the  maximum  AT  which  can  be  reached  is  60°C.  This  is  a 
fundamental  limitation  which  can  be  overcome  in  2  simple  ways: 

-  use  2  stacks  and  corresponding  heat  exchangers  inside  a  double  Vi  wave 
resonator  and  add  serially  the  AT  of  each.  We  have  shown  that  this  approach 
works.  It  has  potential  for  future  developments  of  thermoacoustic 
refrigeration. 

-  increase  the  stack  length  Ax  by  using  a  fluid  where  the  speed  of  sound  is 
higher  than  in  air.  With  air  at  1  atmosphere  we  were  near  the  limit  as  far  as 
AT  max  is  concerned. 

(ii)  Heat  transfer  at  stack-heat  exchanger  interface. 

The  gradual  transport  of  heat  along  the  stack  during  refrigeration  operation  ends 
when  the  symmetry  is  broken  at  each  end  and  hence  a  heat  exchanger  is  needed  at  each  end  to  dispose 
of  the  heat  or  absorb  it.  At  the  cold  end  the  interface  has  to  transfer  heat  Qc  while  at  the  hot  end  the 
heat  transferred  there  is  Qc  +  W,  where  W  is  the  work  done  on  the  system  by  sound.  Since  at  the 
interface  of  stack-heat  exchanger  heat  is  transferred  by  thermal  contact  of  the  cotton  wool  fibers  to 
the  heat  exchangers,  the  contact  thermal  resistance  will  limit  the  flow  of  heat.  A  contact  thermal 
resistance  can  be  defined  as: 


Rco  =  1/hco  A  (7) 

where  hco  =1.25  ks  (m/a)  (P/H)  (8) 

with  ks  being  a  harmonic  mean  thermal  conductivity  for  the  2  solids  in  contact,  o  is  a  measure  of 
surface  roughness  of  the  2  solids,  m  is  related  to  angles  of  contact,  P  is  the  contact  pressure  and  H 
is  the  microhardness  of  the  softer  solid.  For  a  transistor  casing  and  a  nylon  washer  this  resistance  is 
2°C/W  while  for  transistor  in  contact  with  air  it  is  5°C/W.  For  our  cotton  wool  to  heat  exchanger 
interface,  the  thermal  resistance  is  estimated  to  be  Rc0  =  3.5-7°C/W.  For  a  total  heat  flow  of  2  watts 
the  interfaces  can  easily  develop  a  AT  of  7-15°C.  Moreover  closer  examination  of  our  random  stack 
showed  that  it  was  formed  from  3  layers  of  cotton  wool  pressed  together  with  the  fiber  alignment 
perpendicular  to  the  axis  of  heat  transport.  This  explains  why  a  AT  of  25  °C  above  the  heat 
exchanger  anchored  to  room  temperature  was  observed  at  the  hot  end  of  the  stack. 

(iii)  Random  Fiber  Stack. 

An  important  function  of  the  stack  is  the  storage  of  heat  as  it  is  being  shuffled 
from  one  end  of  the  stack  to  the  other.  This  requires  a  large  surface  area;  cotton  wool  is 


exceptionally  well-suited  for  this  task.  Our  stack  has  an  enormous  surface  area,  7450  cm2.  It 
occupies  60%  of  the  stack  volume  the  rest  being  air.  Actually  such  filling  factor  is  too  large  and  it 
should  be  reduced  by  a  factor  of  ~2  to  accommodate  for  the  thermal  penetration  depth  around  each 
fiber. 

Because  our  stack  is  so  short,  Ax  =  0.5  cm,  thermal  conduction  losses  across  it 
are  a  limiting  factor  for  large  AT.  The  total  heat  leak  across  air  and  the  cotton  wool  fibers  (thermal 
conductivity  —  0.04  W/m/  C)  in  the  stack  is  ~0.32  W  for  a  AT  =  50°C;  this  is  a  significant  fraction 
of  the  estimated  1.2  W  cooling  power  of  our  refrigerator. 

For  short  stacks,  much  as  in  our  case,  a  random  fiber  approach  provides  improved 
performance  over  the  fishing  line  spacers  with  mylar  sheets  and  leads  to  simplicity  in  the  construction 
of  the  stack. 


(iv)  Working  Fluid. 

The  reason  for  choosing  air  at  1  atmosphere  as  working  fluid  was  for  simplicity 
and  to  demonstrate  proof  of  concept.  Since  the  speed  of  sound  in  air  is  low  it  determined  the  length 
Ax  of  the  stack  which  limited  the  maximum  AT  that  could  be  achieved.  A  Prandtl  number  of  0.7  led 
to  large  viscous  losses  in  the  boundary  layer  next  to  the  stack  elements.  However  the  low  thermal 
conductivity  of  air,  0.26  mW/cm/°C,  was  useful  for  maintaining  a  large  AT  across  the  stack. 

The  higher  speed  of  sound  in  the  He-Ar  mixture  will  dictate  a  longer  stack  (Ax 
—  0.72  cm)  which  will  allow  a  maximum  AT  of  ~86°C  to  90°C.  Such  improvement  will  have  to  be 
weighted  against  the  heat  losses  by  the  gas  across  the  stack,  being  larger  since  the  mixture  thermal 
conductivity  is  0.6  mW/cm  °C 

(v)  Multiple  Stacks. 

This  approach,  demonstrated  above,  is  important  for  at  least  2  applications: 

the  limiting  AT  achieved  across  a  stack  is  determined  by  the  critical 
gradient  and  the  length  of  the  stack.  By  adding  stacks  in  series 
thermally  a  larger  overall  AT  could  be  achieved,  opening  the  way  for 
very  low  temperature  refrigeration  using  thermoacoustics, 
operation  at  high  frequencies  causes  all  the  dimensions,  including  the 
stack,  to  be  reduced.  Ours  is  0.5  cm  for  5,000  Hz  compared  to  the 
one  of  Wheatley  et  al,  [2]  and  of  Hofler  [5]  which  is  ~7  cm  long. 
Increasing  the  frequency  above  5,000  Hz  would  require  even  shorter 
stack  lengths  Ax  which  is  impractical.  However  multiple  stacks  in 
series  would  overcome  this  problem,  thus  making  it  possible  to  go  to 
the  ultrasonic  range. 

(vi)  High-frequency  operation. 

There  are  many  reasons  why  high  frequency  operation  of  a  thermoacoustic 
refrigerator  can  lead  to  important  developments  of  this  field.  Gifford  and  Longsworth  in  their  classic 
paper  [10]  on  Surface  Heat  Pumping  point  out  that  there  is  an  optimum  pulse  rate  for  heat  pumping 
which  scales  with  the  difiusivity  of  the  gas  and  geometric  factors.  Swift  [7]  in  his  paper  on  Similitude 
in  Thermoacoustics  shows  how  to  scale  in  size  and  frequency  various  systems.  For  example  he 
showns  how  cutting  the  dimensions  of  a  heat  engine  in  half  and  doubling  its  mean  pressure  will  leave 
temperature  and  pressure  ratios  unchanged,  but  cut  heat  and  acoustic  powers  in  half,  and  double  the 


pump  frequency. 

A  more  practical  reason  for  going  to  high  frequencies  is  that  the  coupling  of  a 
piezoelectric  driver  to  air  is  more  efficient  at  high  frequencies  since  the  acoustic  impedance  of  air  is 
higher;  thus  it  provides  a  better  match  to  the  piezo.  For  audio  a  transformer  is  used  in  the  form  of 
2  piezos  face  to  face,  i.e.  the  bimorph  with  parallel  connection  of  its  plates  to  reduce  its  impedance. 
Further  improvements  were  achieved  by  cone  loading  the  piezo.  However,  this  cone  will  not  be 
necessary  when  the  pressure  of  the  working  fluid  is  raised.  This  brings  out  another  advantage  of  high 
frequency  operation  (and  smaller  device):  very  high  fluid  pressures  can  be  used  before  limitations  of 
strength  of  materials  come  into  effect  since  the  surface  area  of  our  device  is  quite  small  (~  90  cm2). 

An  important  consideration  for  high  frequency  operation  of  this  refrigerator  is  that 
large  critical  gradients  VT^  can  be  attained.  Since  this  parameter  is  essentially  T,/xl5  the  temperature 
change  Tx  due  to  the  acoustic  pressure  variation  Pj  and  the  displacement  in  the  sound  wave;  this 
leads  to  a  large  temperature  change  T,  with  small  displacement  Xj  since  xx  -  u/co  (where  Uj  is  the 
particle  speed  in  the  sound  field).  Compression  and  expansion  in  a  sound  field  causes  a  gas 
temperature  oscillation  which  leads  to  a  temperature  difference  between  the  gas  and  the  stack.  Such 
temperature  difference  causes  a  heat  flow: 

Q  =  ^PTmP,x,(o  (7) 


where  X  relates  to  heat  exchange.  On  the  other  hand,  a  temperature  gradient  along  the  stack  causes 
a  reverse  heat  flow  from  stack  to  gas: 


C  pm  V  T  xf  0) 

p  ■  m  ml 


where  pm  is  the  density  of  the  gas,  Cp  its  specific  heat,  Tm  is  the  ambient  temperature,  and  VTm  the 
temperature  gradient  [10].  This  shows  how  a  small  Xj  and  large  lead  to  a  large  temperature 
difference  across  the  stack  and  hence  to  a  low  minimal  temperature. 

High  frequency  operation  also  favors  a  high  power  density.  The  energy  flux  H 
per  unit  volume  is  proportional  to  the  pump  frequency  f  being  [7]: 

H/V  ~(f72)Tm  pP12/(]+ej)(pma2)  (9) 


where  es  is  a  parameter  which  compares  the  thermal  properties  of  the  gas  (pmCpfiJ  to  those  of  the 
stack  (psCp6s),  and  a  is  the  speed  of  sound.  Power  densities  of  ~10  W/cm3  can  be  achieved  at  f 
-5,000  Hz. 

Finally,  high  frequency  operation  for  a  resonant  system  leads  to  small  total  volume 
for  the  refrigerator.  This  is  attractive  from  the  point  of  view  of  applications:  compactness  and  rapid 
cool-down  are  important  factors. 


Having  presented  the  goals,  the  approach,  and  the  experiments  performed  on  a  high 
frequency  thermoacoustic  refrigerator,  it  is  important  to  put  into  evidence  the  quantitative  aspects 
of  this  work  and  all  the  new  things  that  have  emerged  from  this  research.  This  research  relies  on 
theoretical  foundations  laid  down  by  Rott  et  al  [1 1]  and  highlighted  by  Swift  [7], 

(i)  DETAILS. 

In  this  section  some  of  the  details  and  important  qualities  for  the  high 
frequency  refrigerator  will  be  listed  and  compared  to  those  of  a  500  Hz  refrigerator  discussed 
in  ref.  7. 


High  Freq.  Refrigerator 

medium:  air 

f:  5,000  Hz 

overall  length:  10  cm 

stack  length  Ax:  0.5  cm 

normalized  stack  length  Ax/A:  0.02-0.03 

VTcrit:  120  K/crn 

stack  surface  area:  7,550  cm2 

total  perimeter  II:  151  m 

thermal  penetration  depth  6k:  0.04  mm 

acoustic  displacement  amplitude  Xj  (typii 

AT:  50°C 

total  volume:  94  cm3 

sound  pressure  p!  (at  160  dB):  2*1 03  pa 

ambient  pressure  Pm:  1  atmosphere 


Regular  500  Hz  Regrigerator, 

helium 

500  Hz 

~40  cm 

8  cm 

0.04 

15  K/cm 

3920  cm2 

4.9 

0.1  mm 

0.2  mm  2  mm 
90  K 
1200  cm3 

2.9*  104  pa  (at -180  dB) 

10  atmospheres 


For  a  ratio  of  actual  temperature  gradient  to  the  critical  gradient,  r=VTm/VTcrit,  between  0.46  and  0.7 
our  cooling  power  has  been  typically  0.4-0. 8  watt  at  a  sound  level  of  160  dB  at  the  stack.  Slightly 
better  performance  was  achieved  with  a  square  wave  drive  than  with  a  sinesoidal  drive. 

From  the  design  of  the  design  of  the  high  frequency  thermoacoustic 
refrigerator  and  choice  of  components,  high  power  capability  or  demonstrated  by  the  low  frequency 
refrigerator  [12]  was  not  expected.  However  unusual  and  new  features  have  emerged  from  this 
research. 


(ii)  Special  Features  of  High  Frequency  Thermoacoustic  Refrigerator. 

Research  on  the  high  frequency  refrigerator  has  shown  interesting  features 

such  as: 

-  large  critical  temperature  gradients  making  it  possible  to  get  to  very  low 
temperatures. 

-  device  is  mainly  a  pressure  device  rather  than  a  displacement  device.  It 
would  be  an  ideal  device  as  a  second  stage  to  a  regular  low  frequency 
thermoacoustic  refrigerator  making  it  possible  to  achieve  ultralow 
temperatures. 

-  offers  the  possibility  to  use  multiple  stacks  and  hence  multiple  stage 
devices.  Such  devices  could  be  run  in  the  ultrasonic  range. 

-  well  adapted  to  miniaturization. 
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-  simple,  cheap,  compact. 

-  its  small  volume  leads  to  a  short  cool-down  time.  In  fact  from  this  point 
of  view  many  parallel  high  frequency  refrigerators  would  perform  better 
than  a  low  frequency  refrigerator. 

-  attractive,  because  of  compactness  and  simplicity,  for  cooling: 

•  biological  samples 

•  electronic  devices 

•  high  speed  electronics 

•  infra-red  detectors 

-  takes  advantage  of  high  efficiency  of  piezo  drivers. 

-  smaller  heat  dissipation  in  piezo  driver  than  in  electromagnetic  driver. 

-  high  power  density. 

-  large  surface  area  stack. 

-  miniature  heat  exchangers. 

-  well-suited  for  high  pressure  operation. 

(iii)  Further  Developments. 

This  refrigerator  needs  further  developments  in: 

-  heat  transfer  at  stack-heat  exchanger  interface. 

-  application  of  gas  mixtures. 

-  optimum  acoustic  drive  frequency. 

-  limitations  of  refrigerator. 

In  fact,  the  above  areas  are  part  of  the  optimization  steps  discussed  previously 
in  the  development  of  this  refrigerator. 
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-  Lectures:  •  New  Materials  Workshop  in  Luxor,  Egypt,  1997 

“Thermoacoustics” 

•  University  of  J.  Fourier,  Grenoble,  France  1994 
“Cooling  with  Sound” 

•  E.T.H.  Zurich,  Switzerland,  1994 
“Thermoacoustic  Refrigerator” 

•  Brigham  Young  University,  Utah,  1994 
“Thermoacoustic  Refrigerator” 
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